function fP = fP(x,a)
xmin=min(x);
fP=(a-1).*(xmin.^(a-1))./(x.^(a));
end

